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Abstract. - The low energy spectrum of the ferromagnetic Kondo lattice model on a iV-site 
complete graph extended with on-site repulsion is obtained from the underlying spl(2,l) algebra 
properties in the strong coupling limit. The ferromagnetic ground state is realized for 1 and 

J^T , 7V+1 electrons only. We identify the large density of states to be responsible for the suppression 

of the ferromagnetic state and argue that a similar situation is encountered in the Kagome, 

^-! • pyrochlore, and other lattices with flat bands in their one-particle density of states. 

o\ 
o 
o 
o 

It is believed that some aspects of the electronic properties of the strongly correlated 

transition metal oxides, like manganites M , can be revealed by considering the Kondo-lattice 

I . Hamiltonian Hkl =T + Hint , with 

T = -J^Ujclc^ tij>0, (1) 

o 



Hint = Y ^2 C j a S ^ (T a l3 C jf3 + U J2 n 3^ n 3^ ( 2 ) 



using standard notations. The kinetic part T describes the hopping of conduction electrons 



on a lattice (a =1,1 is the spin index). The interaction part Hj n t includes the intra-atomic 
ferromagnetic exchange between the conduction electrons and localized spin S c core electrons M 
with cr a p/2 and S CJ - as their respective spin operators (er denotes a vector of Pauli matrices), 
and the on-site Coulomb repulsion U > between the electrons. While here we are interested 
in Jh > Hund's coupling, the Kondo-lattice Hamiltonian has been extensively studied for 
Jh < 0, as describing heavy fermion systems Q. 

Because the Hund's coupling is typically larger than the hopping, the energetically unfavored 
low spin states are neglected (that is when electron and core spin are antiparallel) , and 

(*)on leave from Res. Inst, for Solid State Physics, Budapest, Hungary. 
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we get the quantum double exchange [|4| M. The next usually used approximation neglects 
the quantum spin fluctuations of the S c core spin described by classical variables (spherical 
angles 9 and <ft). These two approximations lead to the double-exchange model || g, with 
Hamiltonian 7Yde — — X) ^ijflfj > which describes charges as noninteracting spinless fermions 
moving in a disordered background of classical spins with effective hopping amplitudes 



tij tij 



cos — cos -£ + sin — sin -L e K<l>i-<l>i) 
2 2 2 2 



(3) 



The charges can freely propagate provided the core spins are aligned and therefore ferromag- 
netism is favored. The main effect of finite Jh is to introduce antiferromagnetic exchange 
between the core spins, which will hinder the free propagation of the charges, resulting in 
a competition between ferromagnetic and antiferromagnetic ordering. To check the scenario 
presented above, different numerical methods are applied and it has been found that on the 
cubic lattice the ferromagnetism is realized for a wide range of electron concentration pi. 

While the scenario described above is widely accepted, in this paper we would like to point 
out that the structure of the underlying lattice is important and we can expect suppression 
of the ferromagnetic phase when the noninteracting one-particle spectrum has nondispersing 
states forming the so called flat bands, like in Kagome or pyrochlore lattice. This might 
appear surprising as, according to the Stoner criterion, the ferromagnetism goes hand in 
hand with large density of states. As it will turn out at the end, the reason behind this 
phenomenon is rather simple and will be demonstrated as follows: We first derive the low 
energy effective Hamiltonian in the strong coupling limit. Then the underlying spl(2,l) 
dynamical supersymmetry allows us to solve the effective model on the iV-site complete graph 
referred as G(N) and defined by a infinitely long ranged hopping tij = 4(1 — <5y) (on a lattice 
with N sites, each site has N — 1 neighbors). The algebraic approach turns out to be a very 
useful tool, and in particular allow us to show that all the wave functions can be obtained from 
the noninteracting ones by an extended Gutzwiller projection. The breakdown of the standard 
scenario is already present in our toy model: the ground state is singlet apart from the case 
of 1 or N + 1 electrons in the system, and the reason lies in the large density of states in the 
one-particle spectrum of Q(N). Finally, we will present arguments that the same mechanism 
will destabilize the ferromagnetic phase in the lattices with flat bands. 

To derive the effective Hamiltonian we start at t = 0, where the different lattice sites 
decouple. An empty site has energy 0; one electron can form with the core spin either the high 
(S c + 1/2) or low (S c — 1/2) spin state, with energies —JhS c /2 and Jh{S c + l)/2, respectively; 
finally two electrons on a site results in a state with energy U and spin S c . This can be 
summarized by representing the high and low spin states using auxiliary fermions / and d, 
respectively, and the spins by Schwinger bosons b a , so that Sj — J2 a a "ja( <T a.p/'^)^ > i6 an( l 



V2S C + 1 ' 3i V25 c + 1 






The anticommutation relation for electrons requires the constraint ^ Q n b , a — n- + n cl : = 2^ 
to be satisfied at each site and the interaction part becomes diagonal: 

Wtot = 3 -f Y, [( S c + !)»i - S - n j - n j n l] + U Y, n j n r ( 5 ) 

3 3 

Choosing U = Jh/2 we can eliminate the four fermion term n^ni. Now, using standard 
techniques p[ , we can apply a canonical transformation to get the effective Hamiltonian which 
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is the expansion in t/ Jh around the atomic limit and in the lowest energy subspace, where we 
keep the / fermions only, it reads (see also p0[ ) 

JnN e S c sr-^ tij + + 

H cff = ^——J i b la b ] J 3 

i,j,a 

f . f +3 

~ J2 J g (25 c + 1)3 ( 5 ^ Sa P- ff »» <To! P)fi b icfiU b kv b j0fj +°(jr)- ( 6 ) 

The first order term oc £ in the expansion is equivalent to the quantum double-exchange 
Hamiltonian M, H . It competes with the next order term oc t 2 /Jh which is essentially an 
antiferromagnetic interaction between the spins, and there is no need to include higher order 
terms. The formula above is equally valid for S c — 0, where it describes the large-C/ Hubbard 
model, implying that some of the results below apply also to the t-J model. The procedure 
can be repeated for general U, resulting in a more complicated t 2 / Jh term H 0. 

When specializing to complete graph Q{N) the effective Hamiltonian can be written as 



2/ 



2 



Jh{2S c 



' r. ft H 



afi 



where e = t — JhS c /2, and we introduced the operators S = ^ Sj, Y = (S c + 1)N — N e /2, 

F^ = y*- bj a fj and F a = Y]- bj a fj. The nonvanishing anticommutation relations between the 
fermionic operators are 

{F a ,i^}=y<W + So- a £. (8) 

The bosonic spin operator S satisfy the usual su(2) spin algebra and commutes with Y. The 
system closes with the commutation relations between the fermionic and bosonic operators 

[^S] = -E^- [F a ,Y] = -\F a , (9) 



with their conjugate. The set of relations (0-0) define a spl(2,l) graded algebra spanned by S, 
Y, F^ and F a , F^ and F a , while S and Y span respectively the su(2) and u(l) subalgebras p2| . 
Both the electron number and total spin are conserved, but 7i e ff does not commute with the 
F'b - spl(2,l) is not a symmetry of eq. (0). Nevertheless 7i e ff can be expressed in terms of the 
Casimir operators of spl(2,l), su(2) and u(l) - the so called dynamical supersymmetry Once 
the representations [Y, S] of spl(2,l) are decomposed into multiplcts (S, Y) of su(2) x u(l), we 
know|12]] on each multiplet the values of the different Casimir operators of spl(2,l), su(2) and 
u(l) and similarly to the t-model |l^, 14 the spl(2,l) representations give the eigenvalues of 



the effective Hamiltonian. 

The generic irreducible representation (irrep) [Y, S) of the spl(2,l) algebra is 8S dimensional 
and can be labeled by two linearly independent Casimir operators [12]. It contains the (S, Y), 
(S - 1/2, Y + 1/2), (S - 1/2, Y - 1/2) and (5 - 1, Y) spin multiplets. Special cases concern 
the irrep [Y — S,S] which contains only the (S, Y) and (S — 1/2, Y + 1/2) spin multiplets 
(dimension AS + 1) and the irrep [Y, S = 1/2] which do not contain the (S — 1, Y) multiplet. 
Applying operators F a and F^ we can walk between the spin multiplets within an irrep. 
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The (S- 1/2, Y- 1/2), (S,Y), (S-1,Y), and (S-l/2, Y+l/2) multiplets of the irrep [Y, S] 
are eigenstates of ^ Q F^F a with the eigenvalues 2Y — 1, Y + S, Y — S, and 0, respectively. 
As the second order of Ti c ff is also quadratic in F operators and can be expressed in terms of 
the Casimir operators, the spin multiplets (S, Y) in the irrep [Y', S'] are eigenstates of eq. (Q) 
with (increasing) energies: 

/ 2t 2 Q max 4. 1 \ / <?n 

2t 2 gmzx _ g \ / £ma X _ £\ f 3 



2^(5,7) = ^-^-^—-j (*___) .^o^), (11) 

%_j,S4.J](5,y) = A e£ , (13) 

where S™ 21 * = NS c +N e /2 is the maximum total spin. The multiplicity M^ii of the irrep [Y, S] 

in an iV-site system is given by the number of times this irrep is contained in [S c +l/2, S C +1/2] N 
(one site irreps reduce to [S c + 1/2, S c + 1/2]). This is determined from the branching rule]12||, 
leading to the following recursion relation: 

S-\S-v-S c \ 
M (N+1) V- V- M (N) 

IV1 [Y,S+l/2] - 2^1 2-~i M [Y+y-S c -l,S-s+l/2Y 

^={0,1/2} s=-S c ~y 

with Mfy_ N , s ,1/21 si = $S N(S c +i/2) as boundary condition. This formula along with the 
energy definition in eqs. (|l0|)-(|l3|) allows an iterative procedure to build the energy density 
distribution of the model. 

In fig. [j] the energies given by eqs . (|l0[) - (|l3|) are compared with those of the Kondo lattice 
computed by exact diagonalization on a small (N = 4) size. At Jh = oo, the observed lowest 
part of the spectrum is exactly the one predicted from analysis of spl(2,l) representations. Up 
to t/ Jh = 0.05 the effective Hamiltonian energies agree very well with the exact ones. For 
larger t/ Jh values, higher correction terms need to be introduced in order to get a quantitative 
agreement. As the multiplicities of the levels shown at the right of fig. 1 do not depend of t/ Jh 
and are those of the dynamical supersymmetry, these correction can be, in principle, calculable. 
To obtain this solution it is essential that the effective Hamiltonian can be expressed using the 
operators of the spl(2,l) superalgebra which, for finite Jn/t, is possible for U = Jh/2 only 
(see also 



15 



For N e = 1 the (S, Y) spin multiplet of the [Y + 1/2, S + 1/2] is missing and the ground 
state is the highest spin state in the [Y, S] irrep. For 2 < N e < N the t 2 j Jh correction in eq. 
( [l0| ) makes the lowest energy state to be the singlet (S integer) or doublet (S half-integer), 
and the low energy spectrum behaves as S(S + 1), like in the infinite range, antiferromagnetic 
Heisenberg model. The particle-hole transformation can be used for N e > N, and we get 
that for N e = N + 1 the ground state is again the highest spin state (like in the Hubbard 
model fllq, |l7|). To summarize, the model is ferromagnetic for 1 and N + 1 electrons only. 

While the algebraic approach gives the spectrum, it does not tell how to get the wave 
functions. For the case of the t-model it was shown [|18| by explicit construction that some 
of the wave functions can be obtained by Gutzwiller projecting the free fermion ones with 
Pg = Yii0-~ n 3i n j'\)- Actually, more is true: for the quantum double exchange all the large- Jh 
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Fig. 1. - The low energy spectrum of the ferromagnetic Kondo lattice on a 4 site complete graph <5(4) 
with S c — 1/2 and N e — 2 for U ■=_Jh/'2- The solid straight lines show the energy of the effective 
Hamiltonian as given by eqs. (|lCJ)-(|l3|) to be compared to the dotted lines from exact diagonalization. 
The multiplicity and the quantum numbers of the irrep the state belong to is also shown. At the 
upper right corner of the plot we can see states with low spin d fermions to appear. 



wave functions can be obtained by projecting out the d fermions with P = Jl 7 (l — nf) from 
a suitable set of the wave functions of the non-interacting Hamiltonian. To this end, one has 
to consider the states |$) corresponding to the (S — 1/2, Y + 1/2) of the [Y, S] irrep. They 
turn out to be exact eigenstates of the Kondo lattice Hamiltonian for any value of Jh with 
energy N e e (eq. fll3|)). These states satisfy Fa\&) =0 and thus they are the states of the non 
interacting Hamiltonian which do not contain zero momentum electron and are invariant with 
P: C a \<&) = and P|$) = |$) (here C a = J2j c j<y)- Starting from these |4>) and adding zero 
momentum electrons, we obtain eigenfunctions of the Jh = U = model. Next, using the 
identity PC\ = FJ.P, the projected state remains in the same irrep with wave-functions which 
is eigenstate of the large- Jh model. In other words, the extended Gutzwiller projection (as 
P = Pg for S c — 0) is exact for the model on Q{N). 

At this point, it is instructive to study the classical (S c — > oo) and t 2 / Jh = limit of the 
model. In the language of Schwinger bosons &t. sw y/SZcos(0j/2) and b'-, « y / SZsm(8j/2)e l ^ j 
[ fl9| . This immediately leads to the hopping amplitudes (Q) of the double exchange model 
with the one-particle Hamiltonian Tii p — e — t(\c)(c\ + \s)(s\) where |c) = J2j cos "2"/j|0) an d 
\s) = 5^- sin-j-e 1 ^/ |0). If we choose the z-axis to point in the direction of the total core 
spin, the \c) and \s) are orthogonal and eigenvectors of Tii p with energies 



t 



^ + |) 



:(N- 



s~ c 



(14) 



Furthermore there are N — 2 states with energy e. While in the lowest energy state the energy 
is linearly decreasing with S, and the fermions can freely propagate when the core spins are 
parallel, for \s) the tendency is reversed: energy is higher for larger S. Filling the one-particle 
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levels, the spectrum for 2 < N e < N electron is — Nt + N e e, e c + (N e — l)e, e s + (N e — l)e, and 
N e e. These energies are equal 0(1/S C ) to the energies (|To|)-([k|), respectively, and thus the 
correspondence between the semiclassical and quantum spectra for this model is established. 
The absence of the true ferromagnetic state is due to the cancellation of the contributions 
linear in S for N e > 2. Again, the projection operator establishes a relationship between 
Jjj = and Jh — > oo states: \c) = PCj|0) and \s) = PCjjO), and the role of the core spins is 
to act like an effective magnetic field which Zeeman splits the energy of C\\Q) states. 

Now, once we have seen that the ferromagnetism is suppressed in the model, we can search 
the reason of such a phenomenon. One of the peculiarities of our model is the large degeneracy 
of the one-particle spectrum when all the spin are aligned. Then S is maximum and we have 
one state at E = e — tN and iV — 1 at E = e. Turning over one of the core spins the S/S c 
ratio is reduced by 2, increasing the lowest energy e c by t, decreasing the energy of one of the 
levels in the degenerate manifold with the same amount t (actually e s in eq. |l4|) and leaving 
the other N — 2 states of the degenerate manifold with the same energy e. In other words, we 
can gain antiferromagnetic energy without loosing total kinetic energy. For this to happen, 
the existence of the macroscopically large number of nondispersing states is necessary. This 
feature is also present in models defined on lattices of corner shared complete graphs (known 
as line graps, see e.g. g(|). Let us call Q r (M) the M-site complete graphs centered around 
the point r and first concentrate on the family where in D dimensions Q r {D + 1) are sharing 
their corners, and (anti)periodic boundary condition are assumed. In D = 2 it is represented 
by the Kagome and in D = 3 by the pyrochlore lattice. In momentum space they have two 
dispersing bands and just above them D — 1 flat bands. We start to fill the flat bands above 
the electron density nc = 2/(D + 1), and we can expect the mechanism we outlined above 
to act when we flip some of the core spins. For D > 3 we can go even further following 
Ref. |M: the Hamiltonian for t 2 / Jh — can be written as Ti. = £(N e ) +1Z, where the number 
£(N e ) is -J H S c N e /2 - tt(N - N e )(l + S c )/(1 + 2S C ) while K = t/(2S c + lJ^Vl,, 
is a positive semidefinite operator. The summation is over the centers of the Q r constituting 
the lattice with corresponding fermionic operator F a ^ r . Choosing an initial state |$) such that 
F a ,r\&) = for any r, the \*ff) = Tl r Fl r FJ r \®} is an exact ground state of the model with 
energy £(N e ), as TZ\^) = 0. The number of good |$) states is large (e.g. (25 c + 1)^ if it does 
not contain electrons) and l^) will inherit its degeneracy and spin. Our wave function |\P) can 
be visualized as the product of suitable chosen ground states (eq. [l(]) of each graph Q r . For 
finite t 2 J Jh the antiferromagnetic term will split the \^) manifold and a high spin state will 
certainly not be the ground state for electron densities larger than jiq = 4/(D + 1) (we added 
two electrons for each Q r by constructing l^)). The construction can be repeated for other 
line graphs as well. For example when Q T (2D) form a D-dimensional hypercubic lattice [p0|, 
we have one dispersing bands and D — 1 flat bands and, consequently, the critical densities 
are lower, nc = 1/-D and nq = 2/D. It should be noted here that in the cubic lattice the 
ferromagnetism is suppressed near zero and half fillings due to lack of available carriers || ||, 
and it is not the case here, where charges are available, however in dispersionless states. 

From the arguments above the following speculative phase diagram emerges, (i) For densi- 
ties below nc the usual double exchange mechanism will stabilize the ferromagnetic phase in a 
large parameter range, (ii) For densities between nc and nQ preliminary studies indicate a kind 
of ferrimagnetic state, (iii) Finally, if n > nQ the physics is governed by the antiferromagnetic 
term oc t 2 j Jh in the effective Hamiltonian with a possible spin liquid ground state. 

Finally, let us compare our result with the Stoner model. The Stoner model is a molecular 
field model of itinerant electrons, where the magnetization is gained at the expense of the 
kinetic energy, thus large density of states is favored. While certainly we cannot apply a mean 
field model to strongly correlated system, it turns out that in the Hubbard model singularities 



K. Pcnc and R. Lacaze: spl(2,l) dynamical supcrsymmetry etc. 7 

in the density of states can help the ferromagnetism to develop in some particular cases [|22| . 
Therefore one is tempted to generalize the conclusions drawn from the Stoner model. Here we 
have shown that the flat bands (at least if they are above the broad dispersing bands) inhibits 
ferromagnetism in the Kondo lattice. If we look more carefully, the fact that the driving force 
for finite magnetization is the kinetic energy gained by spin alignment, we do not expect from 
the beginning the Stoner criterion to apply. 

To summarize, we have shown that the strong coupling limit of the Kondo lattice with 
infinite long range hoppings can be solved exacly using the underlying dynamical supersym- 
metry. We learned that on this particular lattice (i) an extended Gutwiller projection becomes 
exact, (ii) the ferromagnetic ground state is not favored. Finally, extending our result to more 
general lattices we arrived at an interesting conjecture that large density of states is against 
ferromagnetism in the double exchange model. 



We would like to acknowledge useful discussions with P. Fazekas, F. Mila, E. Muller- 
Hartmann and H. Shiba. 
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